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Abstract 

In this paper we consider the recurrent equation 

for p ^ 1 with / £ C[0, 1] and Ai = y > given. We give conditions on / that 
guarantee the existence of y^ such that the sequence A p with Ai = y^ tends to a 
finite positive limit as p — > oo. 

1. Introduction 

The following problem arose in the joint papers of the first author and Dong Li (see [LS08aJ 
and [LS08bJ). Let / be a continuous real- valued function on [0, 1]. Define the sequence A p 
for p = 1, 2, . . . by 

Vi^E/(^i) a A + h (i) 

g=l ^ ' 

and set Ai = y ^ 0. We shall occasionally write A p (y) to emphasize the dependence of A p 
on the initial value y. It is clear that A p (cy) = c p A p (y). Therefore if A p (y) — > oo as p — > oo 
and c > 1, then A p (y r ) — > oo as p — > oo where y' = cy. On the other hand if A p (y) — > 
and < c < 1, then A p (y') — > 0. Thus there exist y + and y~ such that A p (y) — > oo for 
y G (y + , oo) with y + as small as possible and A p (y) —> for y G (0, with y~ as large 
as possible. It is a natural question whether y + = y~ = y^ and whether A p (y^) — > const 
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as p — ► oo. It is easy to see that this constant must be ( f(x)dx) , and it is our first 



1 i 



o 



assumption that the last integral is positive. It is enough to consider the case / f(x)dx — 1 

o 

because if f(x) = Kf(x) for a constant K, then A p (y) = K^Kpiy). If the answer to our 
question is affirmative then A p (y^) is called the separating solution of ([I]). 

This problem was considered previously in [Li] and [Sin07j . The analysis in [Li] covered 
the case f(x) = 6x 2 — lOx + 4 needed in |LS08a| . The analysis in |Sin07] was based on 
a different idea but unfortunately had a number of gaps. This paper is a modified and 
corrected version of [Si n07] . 

Before we give the assumptions we impose on /, we remark that f(x) and f(l — x) 
produce identical sequences. Therefore the existence of a separating solution depends only 
on fi(x) = f(x) +f(l — x). Of course establishing existence of a solution for / guarantees its 
existence for g if gi = f\. Given fi(x) one can find f(x) so that /(l) = 0. Thus we assume 
that /(l) = without loss of generality. Now we impose the following conditions on /: 

1. /eC 2 [0,l], 

2. f% is positive on [0, 1] fl Q, 

i 

3. all complex a + 1 satisfying j t° h(t )M = 1 have the property that Re, < 0, 



4. a numerical condition to be explained later. 

Observe that an assumption similar to 2 is necessary as A p will vanish for p sufficiently large if 
f\ vanishes on too large a set (e.g., if = 0); Assumption 2 effectively ensures that A p > 

X 

1 f 

for all p. Finally we introduce functions f2{x) = —(xf(x))' and ^(x) = / tf 2 {t)dt. 



x 



Define a p > for p ^ 1 by the condition A p (a p ) = 1; Assumption 2 above makes this 
possible. The strategy of the proof will be to show that a p a m sufficiently rapidly. Take 
positive constants A and B with B < 1 < A and consider the inequalities 

B ^ a p | Op - Op_i| < A/p 2+s ; (2) 

where p is given and 5 G (0, |) will be chosen later and will depend on f\. 

Theorem (Main Theorem). Let f satisfy assumptions 1-3 above. If for some p$ (depending 
on A, B, and f\) the inequalities (j2J) hold for p ^ p , then they are valid for all p ^ 1. 
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Our proof will be inductive. We shall assume ([2D for p ^ r and prove it for p — r + 1. 
This will imply that the limit lim a p = exists and A p (aoo) will be the desired separating 

p— >oo ' 

solution. 

The rest of the paper is structured as follows. In Section 2 we derive a recurrent equation 
for a p . In Section 3 we solve this equation using the inductive hypothesis. The last Section 
consists of numerical analysis and outlines further research on the problem. 

The first author thanks NSF for the financial support, grant DMS N 0600996. 



2. Recurrent Equation for a 



p 

We shall denote absolute constants by C with superscipts in the course of this calculation. 
We have that 

A p+ i(a p+ i) - A p+1 (a p ) = -(A p+1 (a p ) - A p (a p )). (3) 
Put 7 = ^, p 2 = p - pi, i = ^j. Then 

v 



1 1 

A p+1 (a p ) = - E f(j')A Pl (a p )A P2 +i(a p ) = - E f(l')(KM P ) - l)(A P2+1 (a p ) - 1) + 

pi=i pi=i 

+ ~ E /(7')(A Pl (a p ) - 1) + l - /(VXWap) - 1) - \ E /(V) = 
pi=i pi=i pi=i 

= ; E /(-/)( v(°p) - i)(a*,+iK) AcyjcvK) - 1) - ~ e /(V) 
^pi=i p pi=i p pi=i 

A similar formula can be written for A p (a p ): 



1 

A p (a p ) = -— E /(7)(A Pl (a p ) - l)(A P2 (a p ) - 1) + 
P pi=i 

+ ^E /i(7)(A Pl K) - 1) l —r E 



' pi=l ^ pi=l 
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Subtracting A p (a p ) from A p+1 (a p ) we get 

Ap+iK) - Ap(ap) = -/ ( — 7-7 J (Ap(ap) - l)(A 1 (a p ) - 1) 
V \P + 1/ 

P-i 

(A Pl (a p ) - l)(A P2 (a p ) - 1) 



~y \P P- 1 J 

1 

+ "E /(V)( A piK) - l)(Vn(ap) - A P2 (a p )) 
P pi=i 

+ (ViK)-i)+ 

Pi=l ^ / 

+ ^ f^r) - E ^/(V) - -^7/(7)) = E4 



We estimate J^. It will be shown that if is the main term while the others have a 
smaller order of magnitude. This term produces the recurrent equation that we shall analyze 
in Section [XI 

It is readily seen that if = 4 1} , where |4 X) | < The reasoning is as follows. Rewrite 

the term as 

It is clear how to bound the second term in the first factor. The second factor can be written 

as 

p-i / \ / 1 

a p - a p -i \ fp-1 



k=1 x a P-i J \ k 

whence it is easy to see that it is bounded by const • Bp ^+s ■ The estimate for the fourth term 
follows. 



We go on to 

p-i 



\P v — 1 

pi=i 



1(7) (A pi (a p )-1). 



For the first factor in the sum we get 



p 



where | e p \ Si ^3- . The second factor is more complicated and we first rewrite it as 



pi / _ \ 



(2 p Cl pi ^ > J dp (2 pi ^ I Pi 



k=2 



k ' 
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The last term of this expression is not more than (-§) 2 • Multiplying out gives the fol- 
lowing expression: 



r(5) 
p 



^2 7/2(7') a P - a Pl + e ( 3 ) 



pi- 



^ p - 1 a 



and kfKCW^d) 2 



Now we deal with three relatively simple terms. Let us begin with the seventh one: 

p-2 



p-2 

= -E 



pi=i 



Pi 



P \P + 1/ P — 1 VP + 1 

I p-2 r 

^e 



+ 



p-1 



p(p - 1) 
p + l 

Pip - !) 



pi=i 
1 



./ 



Pi 



Pi 



p+l J (p + l) VP + 1 



Pi 



+ e 



f 



(4) 



Pi 



P+ 1 



~f 



[f(l)+lf(l)]dl + ef\ 



Our assumption that /(l) = implies that the last integral vanishes. Thus, Ip 7 ^ = and 



(5) 1 < gg) 



It is easy to see that = and |/p U ''| ^ 
(2) 

To estimate Ip we rewrite it as 



r(6)| / cw 



^ ' = E ( 

pi=i v 



tTvY) + m 
p(p - 1) 



+ e 



(6) 



pi 

E 

,fc=i 



p- 



"pi 



P2 



EP 



k=i 



a„ — a 



P2 



"PI 



The terms in the brackets are bounded by C^ttt and C^ttj. Thus the estimate for this 

J Bp° Bp° 

term becomes Bp t+2S ■ 

Finally for the third term we need to estimate 

A P2+ i(a p ) — A P2 (a p ). 

It is not difficult to see that |A P2+1 (a p ) — A P2 (a p )| ^ (-§) 3 - Combining this with the 

remaining factors gives the bound C^^^ ("§) 5 f° r ^ ■ We have used the fact that 
f(x) ^ C(l — x) in the last step. 

Now we can put the seven terms together and see that 



Ap+i(flp+i) - A 



P+i 



'„ \ - - sr 7/2(70 s - ggi (?) 

Pl=1 P 1 a pi 
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where \sp\ ^ (-g) 5 . A simple calculation gives the recurrent equation 

p 



pi-- 



7/2(7') Qp-Qpi I g (8) 
1 P - 1 a Pi P 



with|4 8) i^^(#) 5 . 



Our objective in this section is to derive a recurrent equation for b p = p 



2 a v —a v -i 



(4) 



Thus 



we rewrite (0J using b p rather than a p . We take a positive integer Q ^ p and get 

Q p 



-P 



.Pi=i pi=Q+i. 



7/2(70 Qp~ Q Pi , _(9) 



p-1 a. 



pi 



Now the sum from 1 to Q gives a contribution bounded by ^^-^Q 1 5 - For the sum from 
Q + 1 to p we observe that 



n 1 



q=pi+l 

The left hand side can be written as 



P 



-y "'Pi 



"Pi 



V -i + p( 10 ) 

2 ^ C Pi 
<?=pi+l y 

with I ep 1 L ' ) I ^ (4) 2 provided Q is chosen sufficiently large and independent of p. Using 



this fact we simplify our equation to 



, \- 7/2(7 \^ 

bp+i = -p 2^ -—^ 2^ 



7/2(70 \ 

2 p ' 



pi= 



q=pi + l 



9' 



(5) 



with l^p 11 '*! ^ (^) 5 • After changing the order of summation we obtain the equation 



u p+i 



1 p 

f q=2 



P 



+ e 



(12) 



(6) 



with l^p 12 "*! ^ ^pnr (^) 5 • This is the equation we set out to solve; it is effectively a linearized 
version of the original equation. 



3. Analysis of the Recurrent Equation 

It will be more advantageous to have a continuous equation rather than a discrete one. To 
this effect we need to define b(x) that would agree with b p when x = p. First we observe 
that ([HD can be written as 
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with a different constant in the estimate for the error term. Now we can extend b as follows: 
set b(x) = with b(x) = on [0, 1). Then we have 



i p 
jb( P y)f 3 (y)dy = ^b g fs(^j 

n 9—2 



with q' G [q, q + 1]. It is easy to see that this sum differs from the one in the recurrent 
equation by not more than ^ Y^q=i V The new error term e(x) will incorporate this term 
as well as e^ 13 \ It is also clear that we need to add lower order corrections to e{x) to ensure 
that b(x) remains constant for non-integral x. 
The equation to solve is now 

i 

b( x ) = J b(xy)f 3 (y)dy + e(x). (7) 
o 

The error term is \s{x)\ ^ ^-jy- (we are dropping the dependence on A and B for now). 
Proposition 1. Let f 3 be given as before and let 




Then all b(x) satisfying (JJJ) with e(x) as above are (possibly infinite) linear combinations of 
elements of {J a& ^{x a , x a log x, . . . , x a log fc_1 x} U {b £ (x)} where k = k(a) denotes the multi- 
plicity of a, and the special solution b e (x) has the property \b e (x)\ ^ 



cm 

„26 ■ 



Proof. This proof can be carried out in a simpler way using the Mellin Transform, but we 
shall stick to the Fourier Transform as it is more common. To this end we set x = e^, 
y = e -", 5(0 = b(e*), F(rj) = -h(e-' n )e" n , E(£) = e(e*). We also extend / 3 to be zero on 
(1, oo). We get 

oo 

B(0= / B(£-ri)F(Ti)dTi + E(Z). 



Taking Fourier Transform of this equation yields 

E(a) 



B a) 



1 - F(a) 



Of course we only require that these are equal as distributions. Now F(—ia) = J t a f 3 (t)dt, 

o 

so we need to look where it attains the value one. It is precisely on the set zE. To invert B, 
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we shall integrate along a countour that goes around points in iH (one can easily see them 
to be isolated) and stays on the real line otherwise. The integral away from the poles will 
give b e {x) and can be bounded as follows. We know that ^ Ce~ 2S ^ (hence the Fourier 

Transform is analytic in a strip centered at the real axis) and that is meromorphic. 

Thus the decay rate for (^ jzj^ j (0 — M e ^) is the same as that for E(x). Integrals near 
poles evaluate to residues at those poles, up to constants. For a simple pole at a' the residue 
is e^ Q . Residues at higher order poles are obtained in the same way. The result is immediate 
once we return to the original variables. □ 

The next proposition will allow us to better understand the structure of E. 

Proposition 2. With notation as above, the set 

Sn{u6C: Re a > cr } 

is finite for each cr > — 1- 

Proof. Let us look only at the real part; in this calculation cr = fj, + iv. We have 
i i 
J cos(v\ogt)t»f 3 (t)dt = -^ J sin{v\ogt)j t {h{t)t» +1 )dt. 



It is clear that the last expression tends to zero uniformly in ji as v — > oo provided fi > cr > 
-1. □ 

This Proposition allows us to study E more carefully. Since 

t i 
fs(t) = hit) - 1 J xh{x)dx and J f\(t)dt = 2, 



1 1 
we always have G E. Set F^a) = J t a fi{t)dt and F 3 (a) = J t a f 3 (t)dt. Then 



o 



Fi{°) 7 = F M 



for cr 7^ 1. This means that it suffices to look for solutions to Fi(a) = 1. It is easy to see 
that Fi(a) < 1 when Re a > 1 even without Assumption 3. It is also clear that ^(1) ^ 1. 
Therefore Assumption 3 effectively says that there are no solutions to F 3 (a) = 1 in the 
strip ^ Rea ^ 1 with the exception of o = 0. Thus cr = is the solution with the 
largest real part. However, this solution is extraneous to our problem because it implies that 
' C/p 2 and thus 



%>-i 



a 



a p 



8 



Ya. G. Sinai and I. Vinogradov 



This is only possible when C = 0, so this solution does not work in our situation. To this 
end we define £' = £ \ {0}. 

Suppose £' is nonempty and let <j\ = max{Recr: a G £'} < 0; it exists by Proposition 1. 
Then choose 5 so that <j\ < —5 < 0. Then the slowest decaying solution b p behaves at worst 
like p ai and 




This means that is the desired separating solution. If £' is empty, define 5 = | and 
o\ = — |. Then the same result holds. It is clear that A and S remain bounded in either 
case. 

4. Numerical Analysis 

While the Main Theorem is true for / satisfying Assumption 3, numerical calculations suggest 
that it is true even without this assumption. First we show how asymptotics work for a 
function satisfying this assumption, say f(x) = 6x 2 — lOx + 4. From Figure 2 we see that 
pb p — > 0, which is expected as £' = for this function. It is reasonable to infer from Figure 
1 that a separating solution exists and that the starting value is approximately 1.412729. 

Next we consider the function f(x) = 9x 8 . It has £' w {-0.234067±2.11581i}. Figure 4 
shows p~ ai b p and cos(Ima logp) (the smaller graph is the cosine). We see that the consecutive 
extrema of the rescaled b p are at about the same absolute heights. In addition, we note that 
zeros of the two functions alternate. Therefore, it is plausible that y4cos(Imcrlogp + B) 
will coincide with our function for sufficiently large p. Figure 3 suggests that a separating 
solution exists and that ~ 2.95072. 

Finally we look at f(x) = 13a; 12 . It has £' ~ {0.105896 ± 1.97567?}, and our theorem 
does not apply in this situation. Nevertheless numerics show that Oqo exists, and its value is 
approximately 3.688371 (see Figure 5). The sequence b p in Figure 6 doesn't seem to follow 
the asymptotic prescribed by a%. It is unclear how to pick S for such a function since the 
error term e p does not have good decay when 8 < 0. Apparently Assumption 3 is not 
necessary for the Main Theorem to hold, but in this case the structure of solutions to (jSJ) is 
unclear. 
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Figure 1: The sequence a p for f(x) = 6x 2 — Wx + 4. 
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Figure 2: The sequence pb p for f(x) = Gx 2 — Wx + 4. 
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Figure 3: The sequence a p for f(x) = 9x 8 . 
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